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The Stabilized Kuramoto-Sivashinsky Equation: A Useful Model For Secondary 
Instabilities and Related Dynamics of Experimental One-Dimensional Cellular Flows 
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We report numerical simulations of one-dimensional cellular solutions of the stabilized Kuramoto- 
Sivashinsky (SKS) equation. This equation offers a range of generic behavior in pattern-forming 
instabilities of moving interfaces, like a host of secondary instabilities or transition toward disorder. 
We compare some of these collective behaviors to those observed in experiments. In particular 
destabilization scenarios of bifurcated states are studied in a spatially semi-extended situation, 
which is common in realistic patterns, but was barely explored so far. 
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An interface that is driven out-of-equilibrium fre- 
quently develops a patterned structure, characterized by 
a spatially periodic array of identical cells. This pattern 
can in turn show a set of various secondary instabilities, 
until a possible transition to disorder [l[ . Such a scenario 
has been encountered in various experiments, such as di- 
rectional solidification [2|, i3i], directional viscous finger- 
ing 0, [1, [3| , Taylor-Dean flow 0] , a locally heated thin 
layer of liquid [al or an array of falling liquid columns 



CouUet and looss have proposed a generic model, 
based on broken symmetries that predicts ten secondary 
instabilities from a primary static periodic cellular struc- 
ture. This model reproduced successfully many features 
of secondary modes associated to broken symmetries on 
the primary pattern, for instance parity-broken (PB) do- 
mains of drifting cells [3] or vacillating-breathing (VB) 
mode leading to out-of-phase oscillations. However, the 
model was built under assumptions of slow- varying space- 
phase variables and thus remained valid only close to 
secondary thresholds. Gil H, 17 1 has recently built an 
extension of this model, that includes possible phase- 
mismatch between the primary static state and the bifur- 
cated one. Therefore, Gil's model could reproduce some 
far-from-secondary-threshold behaviors, such as oscillat- 
ing patches left behind a propagative domain, black soli- 
tons or spatiotemporal disorder. 

Alternative approaches using partial differential equa- 
tions or cellular automata have been proposed, where the 
possible dynamical modes are not explicitly introduced 
in the model, but rather appear via the unstable dynam- 
ics of cellular solutions. An example is the stabilized 
Kuramoto-Sivashinsky (SKS) equation, which is inves- 
tigated numerically in this paper. In its non-stabilized 
form, this equation was first built to reproduce some gen- 
eral phenomena in falling film on inclined substrates [18| 
or in flame- front instabilities [l^. The SKS equation is 
the following: 
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The term —af represents the damping term. It has 
been shown [20] that this equation was one of the simplest 
to capture ubiquitous features of pattern-forming insta- 
bilities in interfacial growing fronts. This equation is also 
known to exhibit spatio-temporal intermittency, i.e. co- 
existence of laminar domains and turbulent patches for 
a large number of cells i2Ll, while secondary bifurcations 
have been found as well [20|, albeit for a small number 
of cells (up to 3). However, very few studies focussed 
on secondary instabilities for an intermediate number of 
cells (typically a few tens). This is the condition under 
which most experimental interfacial patterns are investi- 
gated, and it is expected that collective behaviors with 
both spatial and temporal significance show up. 

In this paper we numerically investigate tertiary bifur- 
cations, i.e. destabilization scenarios on secondary dy- 
namical states. These secondary states are themselves 
results of the destabilization of a primary static periodic 
structure. With this semi-constrained (or semi-extended) 
geometry, the pattern is large enough to show collective 
behaviors, and small enough to allow for the tracking of 
the motion and the shape of a single cell. We aim to find 
a generic comprehensive model for complex states that 
appear to be due to non-trivial mode coupling or finite- 
size effects: (a) oscillating wakes behind a propagative 
domain, (b) an amplitude hole corresponding to a phase 
jump in an extended oscillatory state, (c) oscillations su- 
perimposed on a state of drifting cells before its rupture. 
We show that these tertiary states, amongst others, can 
be reproduced by the SKS equation, demonstrating its 
relevance for phenomena beyond the threshold of sec- 
ondary instabilities. This paper is divided into two parts: 
we give a short description of the numerical method, fol- 
lowed by a set of results in the form of spatio-temporal 
diagrams. Then we comment on similarities and differ- 
ences with experiments. 
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Results The resolution of eq. ([T]) is carried out by a 
pseudo-spectral method. The space derivatives are cal- 
culated in Fourier space: a multiplication of the vector of 
the Fourier coefficients, by the vector of the correspond- 
ing wave- numbers times the complex imaginary unit i, 
gives the first space derivative. Any n*'' space derivative 
is calculated via the same type of multiplication that is 
repeated n times. The time derivative of / is then evalu- 
ated by a finite difference method. The choice of a small 
time-step (typically around 10"'^) is suitable in order 
to avoid convergence problems. An implicit scheme has 
been used as well, but is not necessary here: the cellular 
solutions are smooth enough for a simple explicit scheme 
to work. We employ periodic boundary conditions and 
the initial conditions are fixed through the number of 
cells Nc, using about 25 mesh points per cell. 

Inspired by the method employed in [20| , we prescribe 
initial conditions as combinations of sinusoidal functions 
in order to trigger secondary instabilities from a pri- 
mary periodic, static pattern. Most of the initials con- 
ditions consist of a single wave-number k plus random 
perturbations (typical magnitude 1/100). Some states 
needed specific initial conditions, like for example the 
PB drifting cells, which are obtained from a combination: 
/o(x) = sin(fca;) -I- asin(2fca; -I- ^o). The phase shift is ar- 
bitrarily chosen at 0.5, and the amplitude a is chosen to 
be equal to 0.5. These two quantities can take values in a 
certain range (0.3 to 2 for the absolute value of (pQ, 0.25 to 
1 for a) without changing the phenomenon qualitatively, 
but they slightly influence the kinetic properties of the 
selected states. As the purpose of this study is to seek for 
destabilization scenarios of secondary bifurcated states, 
we opt for a set of parameters (fc, a and Nc) such that 
the system is close to the boundary of existence of the bi- 
furcated states. To compute the temporal evolution of a 
localized domain of PB cells, we chose the following initial 
condition: foix) = sin (fca;(l + |(1 -|- tanh(a::iim — x))) + 
asm {2kx + 0o)) x i(tanh(a;iini — x) + 1). Then between 
X ~ and xum, the cells are asymmetrical and have a 
larger wavelength, which corresponds to the domain of 
stability of drifting cells [20| . 

Figure [T] gives a cartography of the main states ob- 
tained by varying both a and k. The stability of each 
state is checked by runs of relatively long duration time 
of one hundred (corresponding to 10^ time steps): if the 
state is not broken, it is considered as stable. The sym- 
bols stand for the domain boundaries: when these are 
crossed, the initial state undergoes a transition to an- 
other one (sometimes disordered). One of the main rea- 
son for the break-up is the occurrence of the Eckhaus 
instability which delimitates the domain for static 
cells. The break-up of drifting states coincides well with 
the neutral curve of the mode of wavelength 2k. The 
OSC regime is stable within a much narrower strip than 
the one predicted by the linear stability analysis [20| . At 
large a, both drifting and static states undergo a transi- 
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FIG. 1: Stability diagram of the primary and secondary states 
in the SKS equation. The domains for static cells, oscillat- 
ing cells (OSC) drifting cells (DRIFT) and oscillating-drifting 
cells (OSCD) are bounded respectively by open circles, dark 
circles, open squares and dark squares. In empty domains, a 
given state is unstable. The full line stands for the neutral 
curve of the mode of wave number k, and the dashed line 
stands for the neutral curve of the mode 2k. 



tion to a flat front and the cellular pattern vanishes: it 
corresponds to the neutral curve of the mode k. 

Typical spatiotemporal diagrams are depicted on fig. 
[21 and reproduce complex collective behaviors (a-c) rem- 
iniscent to various experiments of patterned interfaces. 
Time runs vertically from bottom to top. These dia- 
grams are obtained by plotting f{x,t) every n*'' time- 
step, added to a vertical shift proportional to the value 
of time. The case (a) represents a state of drifting cells, 
initially traveling at constant speed, which undergo an os- 
cillatory instability in a second step, to ultimately break- 
up and enter a disordered regime. The case (b) is sim- 
ilar to (a), except that each cell oscillate out-of-phase 
of its nearest neighbors. These oscillations are likely to 
get amplified and lead to spatiotemporal disorder, like in 
case (a). The case (c) shows a local domain of drifting 
cells, propagating opposite of drift, and leaving oscilla- 
tory patches behind its trailing edge. 

A second kind of behavior concerns phase imperfec- 
tions in extended oscillatory regimes, fig. [31 These di- 
agrams are obtained when an oscillatory state develops 
and exhibits a so-called 'optical mode' (by analogy to 
eigenmodes in phonons) i.e. neighbor cells oscillating out 
of phase. While such a state can be almost perfectly ho- 
mogeneous if the number of columns is even, it will nec- 
essarily adapt to a non-homogeneous spatial phase for 
an odd number of columns. These diagrams show that 
the phase imperfection can be twofold: it can either be 
sharply localized (top diagram on fig. [3]) or evenly dis- 
patched on the pattern (bottom diagram). In this sec- 
ond case, the dashed line describes the iso-phase of the 
double-period state. The first case exhibits an amplitude 
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FIG. 2: Examples of complex dynamics beyond a secondary 
instability, (a) A global drift that undergoes in-phase oscilla- 
tions Eis a first stage to the break-up towards spatiotemporal 
disorder (fc=0.47, q:=0.095, A^c^lS). (b) A global drift that 
undergoes a bifurcation to out-of-phase oscillations (fc=0.53, 
a=0.04, A^c=12). (c) Oscillating wake left behind a domain 
of drifting cells (extract) (a= 0.15, iVc=18). 



FIG. 3: Imperfect extended oscillatory, period-doubling state, 
caused by an odd number of cells and periodic boundaries. 
Top: localization of an amplitude hole. Bottom: progressive 
phase-shift. Both obtained with: k—0.64, a=0.1, A'^c=19. 



hole, i.e. the amplitude of oscillations vanishes in the 
vicinity of the phase imperfection, enabling the pattern 
to have a sharp phase jump of tt. 

Comparisons with experiments - Discussion We now 
argue that the features depicted above reproduce desta- 
bilization scenarios observed in several experiments. We 
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FIG. 4: Measurements of drifting velocities in the SKS equa- 
tion (a,b,c) and in the experimental pattern of liquid jets (d). 
(a) Versus 1/k for two values of a. (b) Versus 1/a for k—0.47. 
(c) Versus 1/a for various k (close-up around threshold), (d) 
Drift velocities versus flow-rate, for different wavelengths (sil- 
icon oil of viscosity r;=100 cP), 



start with a presentation of quantitative data on the PB 
drifting cells. Figures |3] show general trends of the drift 
speed Vdi measured on a globally extended drifting state 
in both the SKS equation (a-c) and a pattern of falling 
jets (d) [13, The distance along x is chosen such 

as the length of a cell equals one. Figure |3]-a shows the 
dependency of Vd versus 1/k, for two values of a. Fig- 
ures |3]-b and -c show Vd versus 1/a: the general trend 
is a sharp increase of the speed just above a threshold 
value for a. Let us also recall that the drift velocity of 
PB cells in the pattern of columns (V^) increases as the 
square-root of P (the flow-rate per unit length), see fig. 
|3]-d. If one tries to get a similar relationship for the SKS 
equation, it turns out that identifying the control param- 
eter with 1/a, is the most relevant choice (compare figs. 
|4]-b, [4]-c and[4]-d). In the pattern of columns, a similar 
dependance on k was found: the drift speed increases 
with 1/k for most of the conditions (fig. |3]-a). In other 
experiments, the identification of 1/a with the control 
parameter of the PB bifurcation is also straightforward: 
in the printer's instability, it is the rotation speed of the 
internal cylinder [Qj and in directional solidification it is a 
combination of the wavelength A and the pulling velocity 
V that reads: XV^. 

The dynamics depicted in figures [2] are commonly ob- 
served in experiments. Case (a) is reminiscent to various 
experiments 0, ■ It can either appear under tem- 
poral modulations of the control parameter f7| or simply 
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when the control parameter is increased sufficiently be- 
yond the threshold of drifting cells 0, ■ Case (b) has 
only been observed in directional solidification and re- 
ferred to as T-2A O in 0]. Case (c) is observed in both 
directional viscous fingering 0] and an array of liquid 



columns jllL [13| . It shows the link between the parity- 
breaking and vacillating-breathing modes: in this case 
the relaxation at the rear wall of the propagative dila- 
tion wave (drifting cells are dilated compared to static 
cells) leads to oscillations at the trailing edge of the do- 
main. These oscillations are the local counterparts of the 
ones in fig. [3l Two of these destabilizing behaviors are 
reproduced on figs. [5]-a and b for the pattern of liquid 
columns [ill. [l3l|. 

We now examine the phase imperfections in an oscilla- 
tory state (figs. [3|). The two types of phase imperfections 
on an extended oscillating state have been observed in the 
pattern of columns 13] for an odd number of columns 
(whereas an almost homogeneous oscillatory state is ob- 
served for an even number of columns), as shown on figs. 
[5]-c and d. It should be noted, however, that the SKS 
equation failed to reproduce a sustained propagating do- 
main of PB cells. Such a domain is bound to shrink and 
to vanish after a while, at least in the explored parameter 
range. Thus probably, a state of PB cells remains stable 
only for an extended drifting state. 

In conclusion, the presented results show that the SKS 
equation is able to reproduce a set of complex situations 
that occur for some secondary instabilities of pattern- 
forming experiments. We have chosen a semi-extended 
pattern (a few tens of cells), as this is frequently encoun- 
tered in experiments, and we have taken parameters and 
initial conditions in order to trigger further destabiliza- 
tion of secondary bifurcated states. The fact that such 
realistic behaviors are reproduced by a simple generic 
equation like ([T]), provides an interesting perspective for 
studying tertiary bifurcations. This includes a disordered 
regime with occurrences of phase defects (such defects 
appear on top of diagram fig. Oa), already reported in 
[20| for a small number of cells and also obtained for a 
larger number of cells (a few tens) in our simulations. A 
comprehensive study of such a regime based on statistics 
of defect occurrences, would be of great interest. 
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FIG. 5: Various collective behaviors in the pattern of liq- 
uid columns (viscosity r;=100 cP), reproduced by the SKS 
equation, (a) Oscillations superimposed on a state of global 
drifting cells, leading ultimately to disorder, (b) An oscillat- 
ing wake in the trailing edge of a propagating domain, (c) A 
phase defect localized in an oscillatory state, (d) Progressive 
phase-shift in an oscillatory state. 
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